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^ ' Abstract 

q . A detailed study of the Counter-Rotating Model (CRM) for generic 

finite static axially symmetric thin disks with nonzero radial pressure 
is presented. We find a general constraint over the counter-rotating 
tangential velocities needed to cast the surface energy-momentum ten- 
i [ sor of the disk as the superposition of two counter-rotating perfect 



fluids. We also found expressions for the energy density and pressure 
of the counter-rotating fluids. Then we shown that, in general, there 
is not possible to take the two counter-rotating fluids as circulating 
KS ' along geodesies neither take the two counter-rotating tangential ve- 

locities as equal and opposite. An specific example is studied where 
we obtain some CRM with well defined counter-rotating tangential 
velocities and stable against radial perturbations. The CRM obtained 
are in agree with the strong energy condition, but there are regions of 
the disks with negative energy density, in violation of the weak energy 
condition. 
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1 Introduction 

The study of axially symmetric solutions of Einstein field equations corre- 
sponding to disklike configurations of matter has a long history. These were 
first studied by Bonnor and Sackfield QXj , obtaining pressureless static disks, 
and by Morgan and Morgan, obtaining static disks with and without radial 
pressure 0, [J. In connection with gravitational collapse, disks were first 
studied by Chamorro, Gregory and Stewart [|]]. In the last years, disks mod- 
els with radial tension ||, magnetic fields || and magnetic and electric fields 
|7j have been also studied. Several classes of exact solutions of the Einstein 
field equations corresponding to static and stationary thin disks have been 
obtained by different authors ||- 18], with or without radial pressure. 



In the case of static disks without radial pressure, there are two common 
interpretations. The stability of these models can be explained by either 
assuming the existence of hoop stresses or that the particles on the disk plane 
move under the action of their own gravitational field in such a way that as 
many particles move clockwise as counterclockwise. This last interpretation, 
the "Counter- Rotating Model" (CRM), is frequently made since it can be 
invoked to mimic true rotational effects. Even though this interpretation 
can be seen as a device, there are observational evidence of disks made of 
streams of rotating and counter- rotating matter |19|, |20| . 

Usually has been considered that the CRM can be applied only when we 
do not have radial pressure and the azimuthal stress is positive (pressure). 
These conditions, however, are very restrictive and, in many cases, we have 
disks models that only agree with them in a partial region. Thus, the CRM 
will be valid only as a partial interpretation of the corresponding disks. An- 
other, common, assumption is to take the CRM as representing two fluids 
that circulate in opposite directions with the same tangential velocity. As 
we will show in this paper, this is not the case and, in general, the two fluids 
circulate with different velocities. Furthermore, in some cases may not be 
possible to obtain a CRM if the two tangential velocities are taken as equal 
and opposite. Also, commonly is assumed that the two counter-rotating flu- 
ids must be taken as circulating along geodesies. We also will show that this 
is not necessary and that only can be made if the radial pressure is constant. 

The aim of this paper is a detailed study of the CRM for generic finite 
static axially symmetric thin disks with nonzero radial pressure. In Sec. 2 we 
present a summary of the procedure to obtain these thin disks models and 
obtain the surface energy-momentum tensor of the disk. In the next section, 



Sec. 3, we consider the CRM for the disk. We find a general constraint 
over the counter-rotating tangential velocities needed to cast the surface 
energy-momentum tensor of the disk as the superposition of two counter- 
rotating perfect fluids. We also found expressions for the energy density 
and pressure of the counter- rotating fluids. Then we shown that, in general, 
there is not possible to take the two counter-rotating tangential velocities as 
equal and opposite neither take the two counter-rotating fluids as circulating 
along geodesies. In Sec. 4, we consider an specific example where we ob- 
tain some CRM with well defined counter-rotating tangential velocities and 
stable against radial perturbations. The CRM obtained are in agree with 
the strong energy condition, but there are regions of the disks with negative 
energy density, in violation of the weak energy condition. Finally, in Sec. 5, 
we summarize our main results. 

2 Relativistic Static Thin Disks 

In this section we present, following closely reference ||, a summary of the 
procedure to obtain finite static axially symmetric thin disks with nonzero 
radial pressure. The metric can be written as the line element, 

ds 2 = e- 2 *[1l 2 d V 2 + e 2A (dr 2 + dz 2 )} - e 2 *dt 2 , (1) 

where $, A and TZ are functions of r and z only. The Einstein vacuum 
equations for this metric are equivalent to the system 

Krr + K zz = , (2a) 

(ft$ r ),r + (ft$,*),« = , (2b) 

K z k r + K r k z - 2K® tr ® >z - K rz = , (2c) 

ft r A, r - ft, z K,z - ft($, r 2 - $/) + n zz = , (2d) 

where we assume the existence of the second derivatives of the functions $, 
A and 1Z. 

A general solution of the above system can be obtained by the following 
procedure. Let v = r + iz and T{y) any analytical function of v. Then we 



take 

K{r, z) = Re T[v) , (3a) 

Z(r,z) = ImT(u), (3b) 

${r,z) = *(K,Z), (3c) 

A(r,z) = 11(11, Z) + In \F'(v)\ , (3d) 



where ^(TZ, Z) and 11(72., Z) are solutions of the Weyl equations |21 



(7^)^ + (7^, 2 ),2 = , (4a) 

n n = n(^ t7l 2 - ^ 2 ) , (4b) 

U tZ = 27W, n V, z . (4c) 



Is easy to see that the condition of integrability of the system ( fig) - 

is guaranteed by the equation (pi|). Also we can see that this equation 



is equivalent with the Laplace equation in flat three-dimensional space for 
an axially symmetric function and so \1/ can be taken as a solution for an 
appropriated Newtonian source with axial symmetry. Once a solution \1/ is 



known, II is computed from (flbD - ( f4q) and so we obtain, from (3a) - 
a solution of the field equations (^a|) - (|2H| ) . 

Now if we assume that TZ, $ and A are symmetrical functions of z and 
that the first derivatives of the metric tensor are not continuous on the plane 
z = 0, with discontinuity functions 

"afe — 9ab,z\ z=0 + ' 9ab,z\ z=0 - ^ 9ab,z\ z _ + 5 

the Einstein equations yield an energy-momentum tensor T^ = Q b a S(z), 
where 8(z) is the usual Dirac function with support on the disk and 

Ql = \{b az 5 z b - b zz 5 a b + g a % - g zz b a b + ¥ c (g zz 5 a b - g az 5 z )} 

is the distributional energy-momentum tensor. The "true" surface energy- 
momentum tensor (SEMT) of the disk, S^, can be obtained through the 



relation 

S b a = [T b a ds n = e A -*Q b a , (5) 



where ds n = ^fg~^ dz is the "physical measure" of length in the normal to 
the disk direction. For the metric (0]) we obtain 

So° = 2e*- A |A„-2$, 2 +^|, (6a) 

S{ = 2e*- A A, z , (6b) 

Si = 2e*- A {^}, (6c) 

where all the quantities are evaluated at z = + . 

We can write the metric and the SEMT in the canonical forms 

g ab = ~V a V b + W a W b + X a X b + Y a Y b , (7a) 

S ab = aV a V b + Pv W a W b + Pr X a X b , (7b) 

with an orthonormal tetrad ea 6 = {V b , W b , X b , Y b }, where 

V a = e~ # (1,0,0,0) , (8a) 

W a = e — (0,1,0,0), (8b) 

X a = e *-A(o, 0,1,0) , (8c) 

Y a = e*" A (0, 0,1,0) . (8d) 

The energy density, the azimuthal pressure, and the radial pressure are, 
respectively, 

° = -Sq , V v = Si , p r = Si , (9) 

and 

q = a + p v + p r (10) 

is the effective Newtonian density 



3 The Counter-Rotating Model 

We now consider that the SEMT S ab can be written as the superposition of 
two counter-rotating perfect fluids that circulate in opposite directions; that 
is, based on the two-perfect-fluid model of anisotropic fluids |23) , we assume 
that S ab can be cast as 

gab = gab + gab ^ ( n ) 

where S 1 " 6 and S ^ are, respectively, the SEMT of the prograd and retrograd 
counter-rotating fluids. 

Let be U± = (U±, U±, 0, 0) the velocity vectors of the two counter-rotating 
fluids. In order to do the decomposition (|TT|) we project the velocity vectors 
onto the tetrad ea 6 , using the relations |24]] 

(12) 



(13) 



Ul = e\U b ± 


Tja ___ 


U c ± e d 


With the tetrad (§) we can write 






Ul = 


V a + XJ±W a 




± 


>/i-u| ' 




and thus 







V a = y- -^— — ^- + = + , (14a) 



Vl- 


- U'iU+Ul - y/1 - \J' 2 + \J-.U$ 


u + -u_ 


v/i- 


- \J\UI - y/1 - XJlU* 



w a = y \j + + -ij — ' ( 14b ) 

where U± = U±/U± are the tangential velocities of the fluids with respect to 
the tetrad. 

In terms of the metric h ab = g ab — Y a Y b of the hypersurface z = 0, we can 
write the SEMT as 

S ab = (a+ Pr )V a V b + (p v - Pr )W a W b + p r h ab , (15) 



and so, using (|T4T), we obtain 



s 



ab 



+ 



/(u_ 


,u_)(i 


-u 2 + ) 


J7£L7£ 




(u + - 


U_) 2 




/(U+ 


,U+)(1 


-U 2 _) 


fjaf/6 



(U + - U_) 2 
/(U + , U_)(l - Ui)3(l - TJl)*(UlUl + U*Ul] 



(U+ - u. 



+ Pr/i 



afc 



where 



/(Ui,U 2 ) = (<r + p r )UiU 2 + p„ - p r 



(16) 



Thus, in order to cast the SEMT in the form flTT|), the mixed term must be 
absent and so the counter-rotating tangential velocities must be related by 



/(U+,U_) = 



(17) 



where we assume that |U±| 7^ 1. 

Assuming a given choice for the counter-rotating velocities in agree with 



the above relation, we can write the SEMT as (11) with 



where 



S a ± b =(a ± +p ± )U£U b ± + P± h 



-,11b 



(18) 



tf++ p+ = 


" 1 - u 2 + ■ 

.U--U+. 


{{*+Pr)V-}, 


(19a) 


O- + P- = 


" 1 - u 2 _ ■ 

.U+-U-. 


{{a+p r )\J + } , 


(19b) 


(T+ + CT_ = 


a + p r - 


- Pip, 


(19c) 



P++ P- 



Pr 



(19d) 



Note that the counter-rotating energy densities a± and pressures p± are not 
uniquely defined by the above relations, also for definite values of U±. 



Another quantity related with the counter-rotating motion is the specific 
angular momentum of a particle rotating at a radius r, defined as h± = 
9tptpU±. We can write 

h± = ; 20 



This quantity can be used to analyze the stability of the disks against radial 
perturbations. The condition of stability, 

d(h 2 ) , x 

\-t > , (21) 

ar 

is an extension of Rayleigh criteria of stability of a fluid in rest in a gravita- 
tional field; see, for instance |2^ . 



Now we analyze the possibility of a complete determination of the vectors 
U±. As we can see, the constraint ([17]) do not determines U± uniquely, and so 
there is a freedom in the choice of U±. The simplest, common, possibility is 
to take the two counter-rotating tangential velocities as equal and opposite; 
that is, 

U± = ± U , (22) 

so that fllTD is equivalent to 

P V ~Pr 



XJ' 



a + p r 



(23) 



This choice, commonly considered, leads so to a complete determination of 
the velocity vectors U±, however, this can be made only when the above 
expression is positive definite. If it is not the case, we will have a CRM valid 
only in a portion of the disk. 

Another possibility, also commonly assumed, is to take the two counter- 
rotating fluids as circulating along geodesies. Let be u± = U±/U± the angular 
velocities obtained from the geodesic equation for a test particle, 

9u,rV 2 + goo,r = 0, (24) 

so that 

u>± = ± uj , uo = . (25) 

9ll,r 



As the spacetime is static, the two geodesic fluids circulate with equal and 
opposite velocities and so this is a particular case of the above considered. 

In order to see if the geodesic velocities agree with (|ITD, we need to 
compute /(U+, U_). In terms of lj± we get 



Uh 






w± , 



(26) 



and so, using ( [7a] ) and (fH5|), we can write 

A + (a + p r -p<p)B 



/(U+,U. 



9n,rVo 



(27) 



where 



A — gn, r Soo + 9oo,rSn, 



B 



9oo9n,r + 9oo,r9n- 



Using the Einstein equations (|2aD - (|2df) and the expressions (|6al) - (6c) for 
the SEMT we can show that 



/(U+,U. 



K 



n r - K$ tr 



dp r 

dr 



(28) 



that is, the counter-rotating fluids circulate along geodesies only if the radial 
pressure is constant. In the general case, however, /(U+, U_) ^ for fluids 
circulating along geodesies and so it is not possible to obtain a counter- 
rotating model with them. 

As we can see of the above considerations, for disks built from generic 
static axially symmetric metrics, the counter-rotating velocities are not com- 
pletely determined by the constraint (|17|). Thus, the CRM is in general unde- 
termined since the energy density and isotropic pressure can not be explicitly 
written without a knowledge of the counter-rotating tangential velocities. 



4 A Family of Disks with Some Stable CRM 

We will now consider a simple specific example where we can obtain some 
CRM with well defined counter-rotating velocities. In order to obtain finite 



static disks with nonzero radial pressure we consider, following reference 
a solution of (pal) - (|25|) obtained by taking 



W 



a: 



Vu 2 ^! 



(29) 



where a > 0, so representing a thin disk located at z = 0, < r < 1. Now 
we take a simple solution of Weyl equations ( pEa{ ) - Q4"(| ) given by pj 



V(K,Z) 



U(JZ, Z) 



/i 


In 


2A; 




, 2 




M 


In 


2F 





R+ + R_ - 2k' 




R+ + R_ + 2k_ 


j 


(R+ + R_) 2 - 4A; 2 



q- 



1 and R^ 



K 2 



4R+ R_ 

(Z±fc) 2 . 



(30a) 
(30b) 



where \x > 0, k = -\/ 

From the above expressions, and using (|6a|) - (|6cf), we can compute the 
energy density and azimuthal and radial pressures of the disks. We obtain 

1 + /j, 2 r 2 ' 



a 



Pf 



Pr 



Pr 



2/i 



a 



a 



k 2 r 2 



2^,2 



p r = Po[l + k 2 r 2 ] 



1 + fi z r 
l + k 2 r\ 

2Jl(c 2 ~fc 2 )/2fc 2 



(31a) 

(31b) 
(31c) 



where p = 2a(a — k)^^ > and < r < 1. 

We consider, in the first instance, two cases where we obtain simple ex- 
pressions for the SEMT. Let be a = 0, so that p r = p v = 0, and 



4/ie 



flTT/2 



° ~ (1 _ r 2)(M 2 +l)/2 ■ ( 32 ) 

We have a disk of dust with positive energy density, so in agree with all the 



energy conditions |Z7|. On the other hand, the energy density is singular 
at the edge of the disk. We can also see that the constraint (|17]) leads to 
U± = 0, so this is a "true static disk", in the sense that we can not obtain a 
counter-rotating interpretation for it. This case corresponds to the Bonnor 
and Sackfield disk of reference |l|. For the second, simple, case we take 
H = k > 0, so that p r = p v = 2a{a — k) and 

a = -4(a - k) 2 . (33) 
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The disks so are made of perfect fluids with constant energy density and 
pressure. As a < 0, the disks do not agree with the weak energy condition, 
but q > 0, as the strong energy condition requires. For these disks we also 
have that U± = 0, so that we do not have a CRM. These are also "true static 
disks". 

For any other value of a > and \l ^ k, the radial and azimuthal 
pressures, p r and p v , are everywhere positive and well behaved for all the 
values of /i and a. For the effective Newtonian density we obtain 

e=—> (34) 

so that is positive everywhere on the disk. Thus, the disks are attractive, in 
agree with the strong energy condition. On the other hand, is easy to see 
that a < when r = 1, for any value of a and /i, whereas that a > at 
r = only if fi > a. That is, in general, the energy density a is not positive 
everywhere on the disk. 

In order to study the behavior of the energy density and pressures we 
perform a graphical analysis of them. We show, In Fig. [l|, the energy density, 
a, for disks with different values of fi and a. We first plot a for a disk with 
/i = 1.5 and a = 0.5, 0.6, 0.8, 1.1, 1.6, 2, 2.3, 2.8 and 3.5. Then we plot a for 
a disk with /i = 5.5 and the same values of a. As we can see, for some values 
of \x and a the energy density is negative everywhere on the disks, whereas 
that for other combination of the parameters the energy density is positive 
in the central part of the disks, but negative in the edge. We also study a 
for many other values of /i and a and, in all the cases, we obtain a similar 
behavior. 

In Fig. H] we depict the radial pressure, p r , for disks with different values 
of \x and a. We first plot p r for a disk with \i = 1.5 and a = 0.5, 0.6, 0.8, 
1.1, 1.6, 2, 2.3, 2.8 and 3.5. Then we plot p r for a disk with /i = 5.5 and 
the same values of a. As we said above, p r is positive everywhere on the 
disks for all the values of [i and a. Now, for some values of the parameters 
p r have a maximum at r = and then decrease monotonly. On the other 
hand, with other values for fi and a, p r is an increasing function of r. As 
with the energy density, we also study the behavior of the radial pressure for 
many other values of /i and a and all the cases considered presents the same 
characteristics. 

Now we study the behavior of the azimuthal pressure, p^, in Fig. [| As 
in the previous analysis, we first plot p^ for disks with \i = 1.5 and a = 0.5, 

11 



0.6, 0.8, 1.1, 1.6, 2, 2.3, 2.8 and 3.5. Then we plot p v for disks with /x = 5.5 
and the same values of a. We find a behavior like with the radial pressure, 
with maximum at r = for some values of the parameters and increasing 
functions of r for some other values. We also consider many other values of 
/i and a and, as was the case with p r and a, we obtain similar behavior. 

We now consider the CRM for the above disks. As p r is, in general, 
dependent of r, we can not take the two counter-rotating fluids as circulating 
along geodesies. However, we can test the possibility of obtain a well defined 
CRM with equal and opposite velocities. In order to do this, we can compute, 
from ( |23"D and ([31]), the tangential velocity U 2 and obtain 



U (2/i - a) + n(2k 2 - afi)r 2 ' ^ 

and, in order to have a well behaved CRM, we impose the condition 

< U 2 < 1 . (36) 

Again, is easier to do a graphical analysis, and so we study the above relation 
for a lot of combinations of the parameters // and a. In many of the cases we 
obtain functions with strong change in the slope, with regions where U 2 is 
negative and with U > 1. In order to see the kind of behavior that we have, 
we present a sample plot in Fig. |], where we plot U 2 for disks with /i = 0.5 
and a = 0.5, 0.6, 0.8, 1.1, 1.4, 1.7, 2.2, 2.3 and 2.8. However, in some other 
cases we obtain U 2 as well behaved functions of r, everywhere positive and 
increasing, but always with U 2 < 1. An example of these cases is shown in 
Fig. |5], where we plot U 2 for disks with a = 5 and \i = 5, 5.1, 5.2, 5.3, 5.4, 
5.5, 5.6, 5.7, 5.8 and 5.9. 

We can also compute the specific angular momentum of these CRM and, 
using (|20D and ([H]), obtain 

(2fi -a) + \(2fi - a)k 2 - 2a/i 2 ]r 2 ' l ' 

Like with U 2 , the graphical analysis is better and, again, we consider the 
above relation for many different values of fi and a. We also find, in many 
of the cases considered, strong changes in the slope of h 2 as an indication 
of strong instabilities of the CRM against radial perturbations. Also we 
found regions with negative values of h 2 , showing that the CRM can not be 
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applied for these values of the parameters. We shown an example of the cases 
considered in Fig. f|, where we plot h 2 for disks with \i = 0.5 and a = 0.5, 0.6, 
0.8, 1.1, 1.4, 1.7, 2.2, 2.3 and 2.8. We also obtain h 2 as increasing monotonic 
functions of r for some other values of fi and a, so corresponding to stable 
CRM for the disks. A sample of these cases is shown in Fig. [5], when we plot 
h 2 for disks with a = 5 and fi = 5, 5.1, 5.2, 5.3, 5.4, 5.5, 5.6, 5.7, 5.8 and 5.9. 
These disks are the same with well behaved tangential velocities, everywhere 
positive and increasing functions of r. 

Finally, we can compute a + + cr_ and a± + p± for the above disks. Using 
(|T|), (|23|) and @T|), we obtain 



cr + + CT_ = 2 p r 



2^,2 



/! 1 + fi r 

a: 



<?±+P± = Pr 



1 + k 2 r 2 
2/i + a 1 + [i A r 



(38a) 



2a 1 + k 2 r 2 



2^,2" 

(38b) 



We study the above relations for the values of the parameters n and a that 
leads to well behaved tangential velocities and specific angular momentum. 
As a sample, in Fig. |] we plot o+ + ct_ and a± +p± for a disk with a = 5 and 
/i = 5, 5.1, 5.2, 5.3, 5.4, 5.5, 5.6, 5.7, 5.8 and 5.9. As we can see, for these 
values of the parameters, the total energy density of the CRM, <x+ + cr_, is 
positive only in the central region of the disks but negative in the rest having 
a maximum at r = and then decrease monotonly. On the other hand, 
a± + p± is always positive for these cases. We see that the CRM obtained 
are not only stable but in agree with the strong energy condition, although 
not with the weak energy condition. 

5 Discussion 

We presented a detailed study of the Counter- Rotating Model for generic fi- 
nite static axially symmetric thin disks, with nonzero radial pressure. A gen- 
eral constraint over the counter-rotating tangential velocities was obtained, 
needed to cast the surface energy-momentum tensor of the disk in such a way 
that can be interpreted as the superposition of two counter-rotating perfect 
fluids. The constraint obtained is the generalization of the obtained in Jl8| , 



for disks without radial pressure or heat flow, where we only consider counter- 
rotating fluids circulating along geodesies. Also, we obtain expressions for 
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the energy density and pressure of the counter-rotating fluids in terms of the 
energy density and azimuthal and radial pressures of the disk. 

We shown that, in general, there is not possible to take the two counter- 
rotating tangential velocities as equal and opposite neither take the two 
counter-rotating fluids as circulating along geodesies. Thus, for disks built 
from generic static axially symmetric metrics, the counter-rotating velocities 
are not completely determined; that is, the CRM is in general undetermined 
since the energy density and isotropic pressure can not be explicitly written 
without a knowledge of the counter-rotating tangential velocities. 

An specific example was considered of a family of disks where we obtain 
some stable CRM with well defined counter-rotating tangential velocities and 
in agree with the strong energy condition, but with regions of the disks where 
the energy density is negative, so in violation of the weak energy condition. 
We also found some disks of the family presenting strong changes in the slope 
of the tangential velocity and the specific angular momentum, indicating so 
strong instabilities of the CRM, and with negative U 2 and h 2 so that was not 
possible to obtain CRM for these disks. Were found also two cases of "true 
static disks", in the sense of U± = 0, and so there is not a possible CRM 
interpretation. 

The generalization of the Counter- Rotating Model presented here to the 
case of rotating thin disks with or without radial pressure is in consideration. 
Also, the generalization for static and stationary disks with magnetic or 
electric fields is being considered. 
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a : ji — 1.5 



a : \x = 5.5 





Figure 1: We first plot the energy density, a, for a disk with // = 1.5 and 
a = 0.5, 0.6, 0.8, 1.1, 1.6, 2, 2.3, 2.8 and 3.5. Then we plot a for a disk with 
ji = 5.5 and the same values of a. 
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p r : jj, = 5.5 




Figure 2: We first plot the radial pressure, p r , for a disk with /i = 1.5 and 
a = 0.5, 0.6, 0.8, 1.1, 1.6, 2, 2.3, 2.8 and 3.5. Then we plot p r for a disk with 
fi = 5.5 and the same values of a. 
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Figure 3: We first plot the azimuthal pressure, p v , for a disk with \i = 1.5 
and a = 0.5, 0.6, 0.8, 1.1, 1.6, 2, 2.3, 2.8 and 3.5. Then we plot p v for a disk 
with /i = 5.5 and the same values of a. 
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Figure 4: We plot U 2 and h 2 for a disk with [i = 0.5 and a = 0.5, 0.6, 0.8, 
1.1, 1.4, 1.7, 2.2, 2.3 and 2.8. 





Figure 5: We plot U 2 and h 2 for a disk with a = 5 and /j, — 5, 5.1, 5.2, 5.3, 
5.4, 5.5, 5.6, 5.7, 5.8 and 5.9. 




o±+P± 




Figure 6: We plot <j + + o~_ and er± + p± for a disk with a = 5 and \x — 5, 
5.1, 5.2, 5.3, 5.4, 5.5, 5.6, 5.7, 5.8 and 5.9. 
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